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Abstract. 

Zero resistance differential states have been observed in two-dimensional electron gases 
(2DEG) subject to a magnetic field and a strong dc current. In a recent work we presented 
a model to describe the nonlinear transport regime of this phenomenon. From the analysis of 
the differential resistivity and the longitudinal voltage we predicted the formation of negative 
differential resistivity states, although these states are known to be unstable. Based on our 
model, we derive an analytical approximated expression for the Voltage-Current characteristics, 
that captures the main elements of the problem. The result allow us to construct an energy 
functional for the system. In the zero temperature limit, the system presents a quantum phase 
transition, with the control parameter given by the magnetic field. It is noted that above a 
threshold value {B > Bth), the symmetry is spontaneously broken. At sufficiently high magnetic 
field and low temperature the model predicts a phase with a non- vanishing permanent current; 
this is a novel phase that has not been observed so far. 



1. Introduction 

The study of non-equilibrium magnetotransport in higli mobility two-dimensional electron gases 
(2DEG) has acquired great experimental and theoretical interest. Strong magnetoresistance 
oscillations (MIRO) and zero resistance states (ZRS) have been observed in high mobility 
GaAs / AlxGai-xAs heterostructures subject to a magnetic field and to microwave radiation 
[U [2]. Our current understanding of this phenomenon rests upon models that predict the 
existence of negative-resistance states (NRS) yielding an instability that rapidly drives the 
system into a ZRS [3j. Two distinct mechanisms for the generation of NRS are known, one 
is based in the microwave- induced impurity scattering [H El El El El [9] , while the second is linked 
to inelastic processes leading to a non-trivial electron distribution function [101 IHl 112] • 

More recently, Hall field-induced resistance oscillations (HIRO) and zero differential resistance 
states (ZDRS) have been observed in response to strong dc electric current excitations 
[131 [m [15] . The effect of a direct dc current on electron transport can be quite dramatic leading 
to zero differential resistance states (ZDRS) [16]. At low temperature and above a threshold bias 
current the differential resistance vanishes and the longitudinal dc voltage becomes constant. 
Bykov et al.|16j discussed their experimental results following an approach similar to that of 
Andreev et al. [3]. In terms of the differential longitudinal resistivity pxx the stability condition 



reads 

Hence, according to the condition in Eq. ([T]) the 2DES is unstable at negative differential 
resistance. The presence of the ZDRS can be attributed to the formation of negative differential 
resistance states (NDRS) that yields an instability that drives the system into a ZDRS. 



2. Model and results 

In a recent work we presented a model [17\ [T8] to describe the nonlinear response to a direct 
dc current applied to a 2DEG in a strong magnetic field. The model incorporates the exact 
dynamics of two-dimensional damped electrons in the presence of arbitrarily strong magnetic 
and dc electric fields, while the effects of randomly distributed impurities are perturbatively 
added. From the analysis of the differential resistivity and the longitudinal voltage we observe the 
formation of negative differential resistivity states. Both the effects of elastic impurity scattering 
as well as those related to inelastic processes play an important role. The theoretical predictions 
correctly reproduce the main experimental features provided that the inelastic scattering rate 
obeys a temperature dependence, consistent with electron-electron interaction effects. 

The more relevant results are presented here. After the current density is worked out, it splits 
into a Drude and an impurity induced contribution: J = + J^"^p. The Drude contribution is 
given by = Ue ev^, where Ue is the electron density and the drift velocity is given as 

er E - Wc-re^ x E 
^d = — 7- — 2"^ • (2) 

Here, is the unit vector normal to the plane of the system, lOc is the cyclotron frequency, m* is 
the effective mass of the charge carrier, and Ttr is the transport scattering time. It is important 
to point out that the Drude term is not added by hand; instead it explicitly appears in our 
formalism because the exact solution of the Schrodinger equation in the presence of magnetic 
an electric fields is obtained from a unitary transformation [TH] that incorporates the solution 
of the classical equation of motion. 

In the presence of elastic scattering an electron exchanges a momentum q = k — k' with 
the impurity scatterers. Hence, the components of the impurity induced density current can be 
expressed as 

jr = T.J'^[f (^.f ) - / )] ii) > (3) 
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where i = x,y; riimp is the impurity density, and / is the Fermi distribution function evaluated 
at the energies £^^q/2 and £/_,,' -q/2 respectively. The Landau levels are tilted by the electric field 

according to Sf^^k = (a* + ^) + ^^k, with 0;^ = k • v^. Notice that the x component of the 
exchanged momentum Qx = kx — k'^ is equivalent to a hopping (shifting of the guiding center) 
in the y direction Ay = I'^Qx', with Is = y/h/eB the magnetic length. The function G^^/ in Eq. 
([3]) is given by 

MM c 

in the previous equation e^j- is the 2D antisymmetric tensor, A^^/ = [cOg + Wc (a* — Z^')]^ ""^c + 
V{q) is the Fourier component of the scattering potential, and F is the Landau level broadening 
factor. In order to take into account the known fact that the width of LLs depends on the 
magnetic field [T9], we consider F^ = j3{2ijJc/i^Ttr) where /3 is a phenomenological parameter that 
takes into account the difference between the transport scattering time Ttr and the quantum 



scattering time Tq. In the case of short range neutral impurities /3 = 1. The matrix elements 
D^^jj are given in terms of the associated Laguerre polynomial, see |18) . 

In recent work is has been reported that the temperature dependence of nonlinear oscillatory 
magnetoresistance in 2DEG subject to a strong dc electric field can be explained if the quantum 
scattering time Tq incorporates a temperature dependence that is attributed to electron-electron 
interactions [20]. The quantum scattering rate is written as 

- = ^ + — , — = , (5) 
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the parameter A has to be experimentally determined, but it is a constant of the order of 
the unity. For short range neutral scatterers the impurity rate can be estimated as l/rjmp = 
Air'^hanimp/m* , where a is related to the scattering length [S]. 

Let us now consider the nonlinear transport regime. In a current controlled scheme: the 
longitudinal density current is fixed to a constant value Jdc and there is no transverse current, 
Jy. This leads to a set of two implicit equations for the density current 

Jy{Ex,Ey) = 0. (6) 

They represent two implicit equation for the unknown Ex and Ey] the equations can be 
solved following a self-consistent iteration [l^. However, it is verified that for the conditions 
that apply in experiments and in the separated LL regions {tOcTtr 3> 1), the solution of the 
previous equations simplify because the following conditions E^ <^ Ey and J*™''^ <^ erieEy/B 
are simultaneously satisfied. Hence, from the first equation in ([6]) it follows that the leading 
contribution to the Hall electric field is given by the classical result Ey = BJ^c/eUe- Whereas 
the second equation in ^ yields a nonlinear Ex — Jdc relation 

Ex = ^^Jdc + — Jr(E^-^^Ey-^) ■ (7) 

The expression for the impurity induced currents 7**"^ is given in Eq. ([3]). In what follows we 
set Jdc = Jx- 

We had previously presented numerical results in very good agreement with the experimental 
ones, based in Eq. ([7]) [l7j, and in the self consistent solution of equations in ([6]) However 
in this work we point out that an approximated analytical Ex — Jx relation can be obtained. 
As we are interested in the separated LL region the dominant effect is provided by the intra- 
Landau transitions. We include the following assumptions: (1) only the LL nearest to the Fermi 
level contribute. (2) The potential is taken as a constant V{q) = 27rfi'^a/m*, corresponding to 
a short-range impurity scattering, and (3) the integrand in Eq. ([3]) has a dominant peak at 
the value Qmax ~ V^Wn^, this allow us to carry out explicitly the integral over the Laguerre 
polynomials. After the angular integral are evaluated, a simply Ex — Jx relation is obtained 
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the parameters in this equation are given as 
m 



(8) 



Ci = — , 62 = nuimpa — - , 63 = . (9) 
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Eq. dH]) depends on the a temperature through the T-dependence of Tty. In the simplest 
approximation, short range impurity scatterers, we can consider that is proportional to r^. 



consequently the T dependence is dictated by Eq. ([5]) . Utilizing this assumption and the simple 
equation in ([5]) we carried out the comparison with the experimental results, observing a very 
good agreement. In particular the transition from positive to zero differential resistance states 
shows the correct dependence on the temperature [21j . 

In this work we shall concentrate in a simple but interesting limit, we take T = and 
consider the analysis of the stability of the NDRS. Therefore the transport scattering rate 
coincides with the impurity scattering rate: l/rtr = l/Timp = 4:Tr'^ Hariimp / m* . In the region 
of positive differential resistivity pxx = > 0, the relation between between and Jx is 
single-valued and an homogeneous uniform electron density throughout the system is stable. 
On the other hand, in the region in which pxx is negative corresponding to a NDRS, the system 
becomes unstable. An approximated scheme (Maxwell construction) replaces the negative slope 
by a horizontal line, i.e. a ZDRS. The real reason for the appearance of a negative slope in 
the Ex — Jx plot is the implicit restraint of uniform electron density throughout the system. 
However in these regions configurations with a non uniform current distribution turn out to be 
the equilibrium configuration of the system. The simplest possible pattern is a domain wall: two 
parts of the sample carry stable density currents Jxi and Jx2 with the same value of Ex- Hence 
we expect a phase transition defined by a critical iso-B line, that passes through a critical point 
that is a point of inflection of this iso--B line; so both and ^^jf vanish at this point. The 
critical values: Ex = Ec, Jx = Jc and B = Be] at which these conditions are satisfied are readily 
determined as 

■^^ = ^^2 ?^3/2' ^^= = ^2 3/2, ^c/'=6 3/2^ • (10) 

with the values of the constants given by: ^ 40.7, .^2 ~ 161.9, and ^3 20.85. Notice 
that the method follows similar steps those used in the analysis of the critical point of the Van 
der Waals equation [22]. 

Using Eqs. ([9]), and ([TO]) we can express Jc, Ec and Be is terms of the parameters of the 
system. In order to get an estimation we consider: rig = 8 x lO^^m"^, /i = WOm^V/s, a = 0.1, 

imp 



nimr, = 4.5 X IQ^^m 2. So we obtain: 



Jc~0.4A/m, Ec^Q.bV/m, B^^^-UT. (11) 

These values correspond to zero temperature, so are not directly comparable with the 
experimental ones, however the order of magnitude show a reasonable agreement [I6j. 
We now define reduced variables 

= — , Er = — , Br = — . (12) 

Jc -DC 



Hence the critical point is localized at Jr = E^ = B^ = I. Using ([8]) and (jl2p . we readily obtain 
the reduced volt age- current characteristic 



Er = 0.25 J, + ^53/2 



J2 

1 + 3.8^ 



-1/2N 



(13) 



which is "universal", in the sense that it docs not explicitly contains any parameter of the 
system. 

Fig. 1 display a series of plots of reduced electric field Er as a function of the reduced 
longitudinal density current Jr ; for various values of the reduced magnetic field. We clearly 
observe a phase transition: below the threshold Br < 1 the slope of the curve is always positive 
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Figure 1. (Color on line). Reduced longitudinal electric field as function of the reduced 
longitudinal current J^, the values of are displayed near to each line. 



and the system is stable. On the other hand if i?^ > 1 and above the threshold bias current 
[Jr > 1) the differential resistivity becomes negative. 

We can define a scalar Lyapunov functional VF[Jr] = / (Prw [Jr], where 

/Er fJr fjE 

Jr [E'^] dE', = j J'^dJ^. (14) 

Using Eq. ([T3D a simple calculations yields: 



w[Jr] = 0.125 J2 + 1.38 5^/2 



l-log<^2 l + Wl + 3.^ 



Br 



1 + 3.^ 



Br 



-1/2- 



(15) 



One can show that VF[Jr] is a Lyapunov functional, i.e. a non-increasing function of time, so 
its minima are stable steady states. To the expression for w[Jr\ in Eq. ()15p one should add a 
field (current) derivative contribution of the form A|V • Ep, where A introduces a domain wall 
thickness scale that produces a positive contribution to W . Fig. 2 display a series of plots of 
w as a function of Jr for various values of the reduced magnetic field. If Br < 1 the functional 
VF[Jr.] has a single minima at = 0, however when Br > 1 new minima appear at J,. 7^ 0. 
Furthermore we observe that for sufficiently strong magnetic field the minima at J,. = is not 
longer the ground state. In general the system is expected to relax to one of the local minimum 
of W, not necessarily the ground state. However in the presence of noise the system is expected 
to escape from the high-lying minima. This suggest that at sufficiently high magnetic field and 
low temperature the model predicts a phase with a permanent non- vanishing current Jm, where 
Jm is the minima of W . This is a novel phase that has not been observed so far, hence it deserves 
further studies. In future work we will use the present results to analyze in detail the formation 
and structure of the domain walls and also to include the finite temperature effects. 
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Figure 2. (Color on line). A plot of w vs Jr, for values of Br = 0.1, 3, andS. 
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